lapply(pre2016list, function1)
[[1]]

[[2]]

[[3]]

[[4]]

[[5]]

[[6]]

[[7]]

This is all well and good. But I found it difficult to remember which zone went where. So I’ve plotted a reference image to go beside the charts. #Pre-2016 Reference Images

arransubsect <- filter(pcs,substr(label,1,4)=="KA27")
lapply(pre2016list, function2)
[[1]]

[[2]]

[[3]]

[[4]]

[[5]]

[[6]]

[[7]]

But thinking about it, I could plot all the years together like this. #Pre-2016 Individual Zones shown on whole island

lapply(pre2016list, function3)
[[1]]

[[2]]

[[3]]

[[4]]

[[5]]

[[6]]

[[7]]

Post-2016

lapply(post2016list, function1)
[[1]]

[[2]]

[[3]]

[[4]]

[[5]]

[[6]]

[[7]]

lapply(post2016list, function2)
[[1]]

[[2]]

[[3]]

[[4]]

[[5]]

[[6]]

[[7]]

lapply(post2016list, function3)
[[1]]

[[2]]

[[3]]

[[4]]

[[5]]

[[6]]

[[7]]

Plot the percentiles as bar charts.

arransimd %>%
ggplot(aes(x=year, y=Percentile)) +
  geom_bar(stat="identity") +
  facet_wrap('DataZone') +
  labs(title = "Arran SIMD Datazones", x = "Year", y = "Percentile") +
  theme(plot.title = element_text(hjust = 0.5))

Splitting the bar charts up.

Ideally now I’d like to annotate the above data to highlight the 2016 plots, and show where the change in DZ occurs. (I.e draw a polygon around S01011171-S01011177). I don’t know how to do that yet, so what I’ll do now is seperate it into 2 plots.

LS0tCnRpdGxlOiAiUiBOb3RlYm9vayIKb3V0cHV0OiBodG1sX25vdGVib29rCi0tLQoKYGBge3J9CmxhcHBseShwcmUyMDE2bGlzdCwgZnVuY3Rpb24xKQpgYGAKClRoaXMgaXMgYWxsIHdlbGwgYW5kIGdvb2QuIEJ1dCBJIGZvdW5kIGl0IGRpZmZpY3VsdCB0byByZW1lbWJlciB3aGljaCB6b25lIHdlbnQgd2hlcmUuIFNvIEkndmUgcGxvdHRlZCBhIHJlZmVyZW5jZSBpbWFnZSB0byBnbyBiZXNpZGUgdGhlIGNoYXJ0cy4KI1ByZS0yMDE2IFJlZmVyZW5jZSBJbWFnZXMKCmBgYHtyfQphcnJhbnN1YnNlY3QgPC0gZmlsdGVyKHBjcyxzdWJzdHIobGFiZWwsMSw0KT09IktBMjciKQoKbGFwcGx5KHByZTIwMTZsaXN0LCBmdW5jdGlvbjIpCmBgYAoKQnV0IHRoaW5raW5nIGFib3V0IGl0LCBJIGNvdWxkIHBsb3QgYWxsIHRoZSB5ZWFycyB0b2dldGhlciBsaWtlIHRoaXMuCiNQcmUtMjAxNiBJbmRpdmlkdWFsIFpvbmVzIHNob3duIG9uIHdob2xlIGlzbGFuZApgYGB7cn0KbGFwcGx5KHByZTIwMTZsaXN0LCBmdW5jdGlvbjMpCmBgYAoKI1Bvc3QtMjAxNgpgYGB7cn0KbGFwcGx5KHBvc3QyMDE2bGlzdCwgZnVuY3Rpb24xKQpgYGAKCmBgYHtyfQpsYXBwbHkocG9zdDIwMTZsaXN0LCBmdW5jdGlvbjIpCmBgYAoKYGBge3J9CmxhcHBseShwb3N0MjAxNmxpc3QsIGZ1bmN0aW9uMykKYGBgCgojUGxvdCB0aGUgcGVyY2VudGlsZXMgYXMgYmFyIGNoYXJ0cy4KYGBge3J9CmFycmFuc2ltZCAlPiUKZ2dwbG90KGFlcyh4PXllYXIsIHk9UGVyY2VudGlsZSkpICsKICBnZW9tX2JhcihzdGF0PSJpZGVudGl0eSIpICsKICBmYWNldF93cmFwKCdEYXRhWm9uZScpICsKICBsYWJzKHRpdGxlID0gIkFycmFuIFNJTUQgRGF0YXpvbmVzIiwgeCA9ICJZZWFyIiwgeSA9ICJQZXJjZW50aWxlIikgKwogIHRoZW1lKHBsb3QudGl0bGUgPSBlbGVtZW50X3RleHQoaGp1c3QgPSAwLjUpKQpgYGAKCiNTcGxpdHRpbmcgdGhlIGJhciBjaGFydHMgdXAuCklkZWFsbHkgbm93IEknZCBsaWtlIHRvIGFubm90YXRlIHRoZSBhYm92ZSBkYXRhIHRvIGhpZ2hsaWdodCB0aGUgMjAxNiBwbG90cywgYW5kIHNob3cgd2hlcmUgdGhlIGNoYW5nZSBpbiBEWiBvY2N1cnMuIChJLmUgZHJhdyBhIHBvbHlnb24gYXJvdW5kIFMwMTAxMTE3MS1TMDEwMTExNzcpLgpJIGRvbid0IGtub3cgaG93IHRvIGRvIHRoYXQgeWV0LCBzbyB3aGF0IEknbGwgZG8gbm93IGlzIHNlcGVyYXRlIGl0IGludG8gMiBwbG90cy4K